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Abstract: Recently a novel large- N reduction has been proposed as a maximally super- 
symmetric regularization of N = 4 super Yang-Mills theory on R x S s in the planar limit. 
This proposal, if it works, will enable us to study the theory non-perturbatively on a com- 
puter, and hence to test the AdS/CFT correspondence analogously to the recent works on 
the DO-brane system. We provide a nontrivial check of this proposal by performing explicit 
calculations in the large- N reduced model, which is nothing but the so-called plane wave 
matrix model, around a particular stable vacuum corresponding to R x S 3 . At finite tem- 
perature and at weak coupling, we reproduce precisely the deconfinement phase transition 
in the M = 4 super Yang-Mills theory on R x S 3 . This phase transition is considered 
to continue to the strongly coupled regime, where it corresponds to the Hawking-Page 
transition on the AdS side. We also perform calculations around other stable vacua, and 
reproduce the phase transition in super Yang-Mills theory on the corresponding curved 
space-times such as R x S s /Z q and R x S 2 . 
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1. Introduction 

The gauge- gravity duality [1] has been one of the most important subjects in string theory 
over the past decade. The most typical example is the so-called AdS/CFT correspondence 
between type IIB superstring theory on AdS^ x S 5 and N = 4 U(iV) super Yang-Mills 
theory (SYM). Even in this case, however, a complete proof of the duality is still missing 1 
partly because it is a strong/weak duality. The region on the string theory side, where 
(semi-)classical treatments of gravity or string theory is valid, is mapped to the strongly 
coupled region in the planar large-iV limit on the gauge theory side. In order to study 
= 4 SYM in the strongly coupled regime from first principles, one needs to have a 
non-perturbative formulation such as the lattice gauge theory. The problem here is that 
the supersymmetry algebra includes translational symmetry, which is necessarily broken by 
the lattice regularization. In order to restore supersymmetry in the continuum limit, one 
generally has to fine-tune parameters in the lattice action. In fact there are considerable 

1 See refs. [2, 3] for some attempts to prove the AdS/CFT correspondence by using the worldsheet 
approach. Also there are remarkable developments based on the integrability [4] . 



- 1 - 



developments in reducing the number of parameters to be fine-tuned 2 , but any lattice 
formulations of N = 4 SYM proposed so far seem to require fine-tuning at least three 
parameters [12]. There are also Monte Carlo studies of the N = 4 SYM based on reduction 
to matrix quantum mechanics of 6 bosonic commuting matrices [13, 14], which confirmed 
the AdS/CFT correspondence for 1/2 BPS operators. 

Here we are aiming at first-principle calculations in J\f = 4 U(iV) SYM respecting 
supersymmetry maximally. Since we are interested in the planar large-iV limit, we may 
well have a chance to use the idea of the large- iV reduction [15]. It asserts that the planar 
large- N limit of gauge theories can be studied by dimensionally reduced models, which can 
be obtained by dimensional reduction. The original idea does not work in general because 
of the spontaneous breaking of the U(1) D symmetry in the reduced model [16], which led 
to various proposals [16, 17, 18, 19, 20, 21, 22]. Since reduced models can be regularized 
by making N finite, one may avoid introducing the lattice structure in space-time [18], 
which causes the breaking of supersymmetry. However, it seemed rather difficult to avoid 
the problem concerning the instability of the U(1) D symmetric vacuum without breaking 
supersymmetry. 

Let us recall here that N = 4 U(7V) SYM becomes conformally invariant if one sets 
all the moduli parameters (represented by the expectation values of adjoint scalars) to 
zero. Hence the theory on R at the conformally invariant point in the moduli space is 
equivalent to the theory on R x S" 3 through conformal mapping. In fact, R x S 3 is obtained 
as the boundary of AdS^, when one uses the global coordinate. Thus N = 4 SYM on 
R x S 3 appears naturally in the context of the AdS/CFT correspondence, for instance, in 
the so-called pp-wave limit [23, 24, 25] and in the bubbling AdS [26, 27, 28]. 

For our purpose it is intriguing to dimensionally reduce the theory by collapsing the 
S 3 of R x S 3 to a point. The one-dimensional gauge theory obtained in this way is nothing 
but the plane wave matrix model (PWMM) [23] as pointed out by ref. [29]. 3 The PWMM 
can be regarded as a mass-deformation of the Matrix theory [31] preserving maximal su- 
persymmetry, where the mass parameter corresponds to the curvature of the S 3 before 
dimensional reduction. 4 The model possesses many classical vacua, all of which preserve 
maximal supersymmetry. 

Recently it has been conjectured [32] that if one picks up a particular classical vacuum 
of the PWMM, which corresponds to a sequence of fuzzy spheres with different radii, one 
can actually retrieve the theory before dimensional reduction in the planar limit. (See refs. 
[33, 34, 35] for earlier discussions.) This conjecture may be viewed as a new type of large- iV 
reduction, which extends the original proposal for the flat space-time to a curved one, and 
at the same time solves the aforementioned problems concerning the vacuum instability. 

2 See refs. [5, 6, 7, 8, 9, 10, 11] for recent works. 

3 In ref. [29] an equivalence between PWMM around the trivial vacuum and the pp-wave limit of TV = 4 
SYM on R x S 3 has been shown at one loop in the pure scalar sector. However, discrepancies are found 
by four- loop calculations [30] . This connection between PWMM and SYM should not be confused with the 
large- N reduction we are going to discuss. 

4 Originally PWMM appeared as a generalization of the Matrix theory to the pp-wave background [23]. 
It is also often referred to as the BMN matrix model in the literature. While we obtain formally the same 
model in the context of large- iV reduction, the interpretation of the model is different. 
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The classical instability is avoided since the PWMM is a massive theory, and the quantum 
instability is avoided, too, since the vacuum preserves maximal supersymmetry. Since the 
planar limit is taken in the reduced model, the instanton transition to other vacua is also 
suppressed and the "fuzziness" of the spheres is removed. 

Viewed as a regularization of the N = 4 SYM on R x S 3 , the present formulation 
respects the maximal SU(2|4) supersymmetry (with 16 supercharges) of the PWMM, and in 
the large- N limit the symmetry is expected to enhance to the full superconformal SU(2, 2|4) 
symmetry, which has 32 supercharges. Considering that the conformal symmetry is broken 
by any kind of UV regularizations, this regularization is optimal from the viewpoint of 
preserving supersymmetries. 

Let us also emphasize that if one naively regularizes the N = 4 SYM on R x S 3 by 
introducing an upper bound on the angular momenta on S 3 , one necessarily breaks the 
gauge symmetry as well as supersymmetry. This problem can be dealt with in perturbative 
calculations by adding appropriate counter-terms [36, 37], but it is not clear what to do 
with it in nonperturbative calculations that we are aiming at eventually. In the present 
formulation, the size of the matrices plays the role of the ultraviolet cutoff, which neatly 
respects both gauge symmetry and supersymmetry. 

In this paper we test the novel large- N reduction at weak coupling. In supersymmet- 
ric theories, it often occurs that certain properties in the strongly coupled regime remain 
qualitatively the same in the weakly coupled regime. For instance, the AdS/CFT corre- 
spondence at finite temperature [38] suggests that there is a first-order phase transition in 
the strongly coupled regime of N = 4 SYM on Rx S 3 in the planar limit, which corresponds, 
on the gravity side, to the Hawking-Page transition [39] between the AdS space-time and 
the AdS black hole. In fact, even in the weak coupling limit of the N = 4 SYM, there 
exists a first-order deconfinement phase transition [40, 41], which is conjectured to be a 
continuation of the one in the strongly coupled regime. 5 We confirm the novel large-iV 
reduction at weak coupling by showing that the PWMM indeed reproduces precisely the 
above phase transition. The main results of this work were reported briefly in our previous 
publication [44]. As a related work, a test of the large- iV reduction has been performed in 
the high temperature limit up to two-loop [45] (See also ref. [46].). An application of the 
large- N reduction to M = 1 SYM on R x S 3 was discussed in ref. [47] . 

This paper is organized as follows. In section 2 we briefly review the deconfinement 
phase transition in J\f = 4 SYM on R x S 3 in the weak coupling limit. In section 3 we de- 
scribe the large- iV reduction proposed in ref. [32]. In section 4 we discuss the weak coupling 
limit of the PWMM around the vacuum corresponding to R x S 3 at finite temperature. 
In section 5 we show analytically that the critical temperature of the PWMM agrees with 
that of N = 4 SYM on Rx S 3 . Derivation of some equations is given in appendix A. 
In section 6 we show analytically that the free energy of the PWMM agrees with that of 
N = 4 SYM on R x S 3 under some assumption. In section 7 we perform Monte Carlo 
simulations to verify this assumption and to demonstrate the agreement of the free energy 
explicitly. In section 8 we show that N = 4 SYM on a more general space-time Rx S 3 /Z q 

5 Thermodynamical properties of TV = 4 SYM on R x S 3 have been studied also in refs. [42, 43] 
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can be obtained by choosing a different classical vacuum of the same model. In appendix B 
we discuss a simpler case of Af = 8 SYM on R x S 2 , which can be obtained from coinciding 
fuzzy spheres. Section 9 is devoted to a summary and discussions. 



2. Brief review of weakly coupled M = 4 SYM on R x S* 3 

In this section we briefly review the calculation [40, 41], which showed that weakly coupled 
J\f = 4 U(/c) SYM on R x S" 3 undergoes a deconfinement phase transition in the planar 
large-fc limit at finite temperature. We will see later that equations similar to the ones 
that appear below are reproduced from the reduced model. In order to make the similarity 
clearer, we use k instead of N for the gauge group in this section. 

Let us introduce a finite temperature T by compactifying the Euclidean time t to 
a circle with the circumference T _1 . Unlike the T = case, the holonomy along the t 
direction becomes nontrivial and it is represented by the holonomy matrix U. One can 
choose a gauge, in which U takes the diagonal form 

[/ = diag(e iai , • • • ,e iak ) , (2.1) 

where a a £ (— tt,tt] (a = 1, ■ ■ ■ , k) are constant in space-time. The a a variables are called 
the gauge field moduli, since they are massless zero modes. 

In the weak coupling limit, all the fields except the gauge field moduli can be integrated 
out at one loop. Since we are going to take the large-/c limit, it is convenient to introduce 
the distribution of the gauge field moduli 

k 

P(0) = lj2S(9-a a ) . (2.2) 

(2=1 

The resulting effective theory for p(9) is given by [40, 41] 

S = k 2 J d6d6' P {6) V(9 - 6') p{6') , (2.3) 

oo 

V{0) = Y,Vp<x*W) , (2-4) 
P =i 

V p = l{l - 6z s (x^) - z v ( X P) - 4(-lf +1 z f { X v)} , (2.5) 
where we have introduced dimensionless parameters 

with R s -6 being the radius of S 3 . We have also introduced the functions 

x + x 2 6x 2 - 2x 3 4xi 

*<(*) = (ir^s. = -(1^)3- . ^ (x) = (T^' (2 - 7) 
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which can be interpreted as the single-particle partition functions for the scalars, the vector 
and the fermions, respectively. Then one can obtain the distribution p{9) exactly in the 
large-/c limit by solving the saddle-point equation 



<W v'{9 - e') P (e') = o . (2.8) 



Obviously the uniform distribution is always a solution to the saddle-point equation. At 
low temperature, it gives the absolute minimum of the effective action. As a consequence, 
the center invariance is unbroken and the free energy normalized by which we call the 
normalized free energy in what follows, vanishes in the large-A; limit. This phase can be 
interpreted as the confined phase. One can show that there is a first order phase transition 
at a critical point determined by 

Vi = , (2.9) 

which gives x c = 7 — 4 a/3 [40, 41] in terms of the dimensionless parameter (2.6). 

Above the critical temperature, the dominant solution has a compact support [—9q, 0q] 
with 6*o < 7r, and the equation (2.8) is satisfied only for 9 E [—9®, 9q]. The center invariance 
is broken and the normalized free energy takes a negative value. This phase can be inter- 
preted as the deconfined phase. Near the critical temperature, in particular, the explicit 
form of p{9) is given by the Gross- Witten form [40, 41] 

p{6) = { ™ ( COS i) V w_sin2 i for W - °° (2 10) 

1 for |6-| > 6-o , 



where 9$ = 2 sin 1 Juj , uj = 1 — \ — — . (2.11) 

V 1 - Vi 

The normalized free energy above the critical temperature is obtained for R S 3 = 1 as [41] 



5 , 



k 2 



0.9877(T - T c ) - 4.248(T - T c )2 - 11.696(T - T c ) 2 + 0((T - T c )5) , (2.12) 



where T c = — l/ln(7 — 4a/3) = 0.37966- • • . This phase transition is speculated to be a 
continuation of the conjectured phase transition at strong coupling, which corresponds to 
the Hawking-Page transition [39] according to the AdS/CFT correspondence [38]. 

3. Large- N reduction for M = 4 SYM on R x S 3 

In order to regularize N = 4 U(A) SYM on R x S 3 respecting supersymmetry maximally, 
we use the idea of the large-A reduction. For that we dimensionally reduce the theory by 
collapsing the S 3 to a point. Thus we obtain a one-dimensional gauge theory 

Spwmm = \ f dttr X -(p t X M f - -AX M , X n ? + ^* f A* - ^h M [X M , *] 
J * ~- 

+^(Xtf + ^{X a ) 2 + ipe l]k XX ] X k + i^ ll2 ^ , (3.1) 
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which is nothing but the PWMM [23]. Here the covariant derivative is defined by Dt = 
dt — i[A, ■ ], where A, as well as Xm and ^, is an N x N matrix depending on t. The 
range of indices is given by 1 < M, N < 9, 1 < i, j, k < 3 and 4 < a < 9. The model has 
SU(2|4) supersymmetry, which includes 16 supercharges, and for fi = it reduces to the 
DO-brane effective theory or the Matrix theory [31]. 

In fact the model possesses many vacua representing multi fuzzy spheres. Explicitly, 
they are given by 



V 

Xi = ^(4 ll) ®l k ) , (3.2) 



7=1 

where are the n-dimensional irreducible representation of the SU(2) generators obeying 
[L^\l^\ = ieijkL^\ The parameters n/ and kj in (3.2) have to satisfy the relation 
Si=i njki = N . All of these vacua preserve the SU(2|4) supersymmetry, and they are 
degenerate. 

In order to retrieve M = 4 SYM on R x S 3 in the planar limit, one has to pick up a 
particular background from (3.2), and consider the theory (3.1) around it. Let us consider 
the case 

kj = k , ni = n + I ^ — for I = 1, ■ ■ ■ , v (3.3) 

with odd v, and take the large- ./V limit in such a way that 

n — ► oo, z^^oo, fc^oo, n >oo 

2 

a 2 k 

with A = - — fixed . (3.4) 
n 

Then the resulting theory is claimed [32] to be equivalent to the planar limit of N = 4 SYM 
on R x S 3 with the radius of S 3 and the 't Hooft coupling constant given, respectively, by 6 

2 

R S 3 = - , X RxS 3 = XV S 3 , (3.5) 
A* 

where V53 = 2tt 2 (R S 3) 3 is the volume of S 3 . 

This equivalence may be viewed as an extension of the Eguchi-Kawai equivalence [15] 
to a curved space-time. It is crucial that we do not need to do anything like momentum 
quenching [16, 18]. As a consequence, the formulation preserves the SU(2|4) supersymmetry 
and the gauge symmetry. For a brief review of the equivalence, see ref. [44]. In the following 
sections, we give a nontrivial test of this proposal by studying the model at weak coupling. 



4. Effective theory for the gauge field moduli 

Let us first consider a perturbative expansion of the PWMM around the most general 
background (3.2). The calculation is analogous to what is done in N = 4 SYM on R x S 3 

6 The relationship between the radius of S 3 and the parameter /i of the PWMM agrees with the one 
obtained in dimensionally reducing M = 4 SYM on R x S [i to arrive at the PWMM (3.1). 
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in section 2. In particular, we introduce finite temperature T in (3.1), and choose a gauge 
in which the holonomy matrix U is diagonal 

V 

=0(ln,®tfj) , (4-1) 
1=1 

Ui =diagfe ta i ,e k i J for 7 = 1, ■ ■ ■ , v , (4.2) 

where a'p G (— 7r, 7r] (a = 1, • • • , fcf). In what follows we use the dimensionless parameters 
(5 and x as defined in eq. (2.6) with R s :t = 2/ ji anticipating the relationship (3.5). 

The effective action for the gauge field moduli around the general background (3.2) is 
obtained in the weak coupling limit [48]. It can be decomposed into 

S'efr = 6S S + S V + 45/ + S v , (4.3) 

where S s , S v , Sf and Sy represent the contribution of a scalar, a vector, a fermion and the 
Vandermonde determinant, respectively. Explicitly, they are given as 

Ss=ir Si^ , S v =± , S f = ± , S V = ± S? , (4.4) 

i,j=i i,j=i i,j=i i=i 



s 



V 



fej kj (n/+nj)/2-l 



^' J) = 2 E E E (21 + l)ln{l + e" W« - 2e"^+ 1 ) cos(a« - <#>)} , (4.5) 

a=l 6=1 J=|n/-nj|/2 

(n/+nj)/2-2 

^ {21 + 1) ln{l + e-^ 2 ^ 2 ) - 2e -^ 2i + 2 ) cos(ai / ) - a[ J) )} 

l=\n I -nj\/2-l+5 IJ 



k, k 



^ j) 4EE 



a=l 6=1 



(n/+nj)/2 

+ ^ (2^ + l)ln{l + e- 2 ' 3 - 2 '-2 e -' 3 - 2 'cos(ai / )-a[ J) )} 

i=|Tij-nj|/2+l 



(4.6) 



fej fcj 



^EE 

a=l 6=1 



(n/+nj)/2-3/2 

£ (21 + 1) ln{l + e~ WD + 2e"^+l) cos^ - ^ J) )} 

l=\m-nj\/2-l/1 
(n/+nj)/2-l/2 

+ £ ( 2Z + X ) ln { X + e~ 2/3(2 ' + ^ + 2e^+5) cos^M - a[ J) )} 

l=\m-nj\/2+l/2 

ki 

E ln 



a^6 



sm ■ 



ai 7) - a|, 7) 



(fe| - fc/)ln2 . 



, (4-7) 
(4.8) 



The second term in (4.8) is needed to make the free energy vanish in the low temperature 
phase. With this constant term, the total effective action (4.3) agrees exactly with the one 
presented in ref. [48]. 

Let us restrict ourselves here to the particular case (3.3), and rewrite the effective 
action (4.3) in a form analogous to (2.3)~(2.5), which is useful for analytical studies. For 
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instance, si 1 '^ in (4.5) can be rewritten as 

k (n/+nj)/2-l 



Si/^ = I £ E + 1) ln{l - e -^(a+i)+*(«i J) -^ J) )} + c.c. 

a, 6=1 (=|nj — nj\/2 

oo fc (n/+nj)/2-l 

= - E 1 E - «i J) )} E ( 2/ + 

p=l ^ a,b=l l=\m-nj\/2 

~ 1 * f fJ) , 9 / e -p^(|n7-nj|+l)( 1 _ e -p^(n7+nj-|nj-n l/ |))\ 

= E^E cos {^ } }^ 1— 7^ • 

p=l y a ,b=l r \ / 



We also introduce the distribution function for ctp as 

k 

P {I \e) = \Y J ^0-a^) for 7=1,... ,u . (4.9) 

a=l 

Then we find that the effective action for the gauge field moduli can be written as 

S cS =k 2 [ dOdO' (0) V {I ' J) (6 - 6') p {J) (6') , (4.10) 
i, j=i ^ 



oo 



V {I > J \6) =^2v^ J) cos( P 6) , (4.11) 



P =i 



V^^ =^{Su - Qz^ J \x p ) - 4 7 ' J )(x p ) - A{-lf +l zf J \x p )} . (4.12) 
We have introduced the functions 

f/ n. , 8 ( X \ ni - nj \ +1 (l -X niJ )\ tt , 

(u) ( \-2 d ( x\ ni - nj \- 1+25lJ {l - x niJ - 25lJ )\ d fx\ ni - nj \ +3 (l-x niJ )\ 

z v • w -x ^ ^—^ + dx~\ r^p ' 



f/n, s 3 9 / x l^- n -'l +2<5 "(l-a; n "- 2,5 ")\ i <9 / x l ri ^ n JI+ 2 (l - x n ") 



(4.14) 



(4.15) 



with rijj = nj + nj — |nj — nj\, which can be interpreted as the single-particle partition 
functions for the scalars, the vector and the fermions, respectively, as in the N = 4 SYM. 
We can analyze the effective theory (4.10) by the saddle-point method in the large-A; limit. 
The free energy is given by 

fpwMM = ~TS eS , (4.16) 
where S e g is evaluated at the dominant saddle point. We will show in section 6 that 

P^^PWMM = p-^SYM (4-17) 
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in the limit (3.4), where the right-hand side is the normalized free energy of the N = 4 
U(fc) SYM on R x S 3 . Therefore, we define the normalized free energy for the PWMM by 
the left-hand side of eq. (4.17). The appearance of the £ factor in the above relation is 
consistent with diagrammatic considerations [32]. 

Since p^\0) = p^\—9) due to symmetry, p^\0) can be expanded as 

1 1 00 

P {I) (e) = - + -Y J P { P^ S (p9). (4.18) 
p=i 

In terms of pp\ the effective action (4.10) is expressed as 

V oo 

^^ 2 EE# , ^ ) # ) . (4-19) 
/,J=lp=l 

In the low temperature regime, Vp 1 '^ given by (4.12) are positive definite matrices. The 
action (4.19) is therefore minimized by a configuration with p p ^ = for all p > 1, which 
implies that ai^ distribute uniformly. The center invariance is unbroken and the normal- 
ized free energy vanishes. This phase can be interpreted as the confined phase. Above the 
temperature determined by 

detV p (I ' J) = for certain p , (4.20) 

we obtain a configuration with pj^ / for some / and p > 1, which implies that some 
of the distributions p^ (6) become non-uniform. As a consequence, the center invariance 
is broken, and the normalized free energy becomes negative. The high temperature phase 
can be interpreted as the deconfined phase. The distributions p^\0) are determined by 
the saddle-point equation 

j d9'V^' J \9 - 9')p (J \9') = for 9 £ [-9^,9^] (4.21) 

j=i 

derived from (4.10), where [— 9q\ 9q^] (9q^ < it) represents the support of p^\9). 



5. Agreement of the critical temperature 

In this section we show analytically that the critical temperature of the PWMM agrees 
with that of N = 4 SYM on R x S 3 . 

Let us first take the n — > oo limit in (3.4). Then, the single-particle partition functions 
(4.13), (4.14) and (4.15) reduce to 
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Note here that zf' J \x) (i = s,v,f) and hence Vp 1 '^ are v x v Toeplitz matrices; i.e., 

z (i+i,J+i) (x) =z (i,J) (x) for i = s, v, / , (5.4) 
pKJ+M+i) = vK/,J) . (5 .5) 

Therefore, we can represent them by zf' J \x) = ^""^(a;), Vp 1 '^ = Vp 1 ^^ and make a 
Fourier transformation as 

oo 

Zi (x,\)= zf\x)e iKX for i = s,v,f, (5.6) 

ftT=-oo 

oo 1 

Vp(A)= £ V^e* KX = -{l-6z s (xP,\)-z v (xV,\)-4(-ir +1 z f (xV,\)} • (5.7) 

Let Up / ' ) be the eigenvalues of Vp I,J \ Then, a theorem for a Toeplitz matrix [49] implies 

lim min-ui^ = minK,(A) (5.8) 

v->oo I p A 

for each p. Since Vp(A) is positive in the low temperature phase, (5.8) allows us to replace 
the condition (4.20) by 

min Vp(A) = for certain p . (5.9) 
A 

This can be further replaced by 

mm{pV p (\)\ = . (5.10) 
A, p I J 

In the appendix A we show that the left-hand side of this equation is given by 

m\pVp(X)\ = Vi(0) = Vi , (5.11) 



mm< 

A, 



where V\ is defined in (2.5). Therefore, the condition (4.20) is indeed equivalent to (2.9), 
which implies that the critical temperature should agree with M = 4 SYM on R x S 3 . 

Let us comment that eqs. (5.1), (5. 2) and (5.3) actually coincide with the single-particle 
partition functions for N = 8 SYM on R x S 2 around a specific monopole background. 
What we have seen in the previous paragraph can therefore be understood as the large- N 
equivalence between N = 4 SYM on R x S 3 and N = 8 SYM on R x S 2 around the 
corresponding multi-monopole background [32]. The agreement of the critical temperature 
for these two theories has been observed numerically in an earlier work [50]. 
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6. Agreement of the free energy 

In this section we demonstrate the agreement of the free energy (4.17). For that purpose 
we pay attention to the single partition functions that appear in the kernel (4.11) of the 
effective action. We find that z\ I,J \x) for the PWMM defined in eqs. (4.13)~(4.15) and 
Zi(x) for the SYM defined in eq. (2.7) are related to each other as 



= z l {x) + Az\ 1] {x) for i = s,v,f, (6.1) 

J=i 

where AzP( x ) = - X JL j_J!L^A<'>(z)} , (6.2) 

A«(x) = Y3^{ x/ + + ^""^(l - ar")} • (6-5) 

For instance, (6.1) for the i = s case can be shown as 

£ ^ ,J) ^ = *#■ (r=^ X> |/_J| - x 2n - v - l+I+J )\ = zs(x) + a4 ; »(x) . (6.6) 

J=l x I x J=l J 

Using (6.1), we find that 

V 

V^ 1 ^ {9) = V{6) + AVW (9) , (6.7) 

j=i 

AVW(6) = Yj -{^ a 4 I) (x p ) + AzP + 4{-l) p+1 Az { f I] (x p )} cos{ P 6) , (6.8) 

where V{9) is the kernel (2.4) for the N = 4 SYM, and the remaining /-dependent part 
AV W (6) decreases exponentially as one moves away from the edges 1 = 1 and I = v. 
Thus we may naturally expect the solution to (4.21) to be 

pW(6)=p(6) + ^P {I) (0) , (6-9) 

where p{0) is the solution for the J\f = 4 SYM, and Ap^\9) decreases exponentially as 
one moves away from the edges. We will see such a behavior explicitly by Monte Carlo 
simulation in section 7. Thus we have 

£l>pW(0)|=0Q • ( 6 - 10 ) 
i=i 

Substituting (6.9) into (4.10) and using (6.10), we obtain 

^ S cS = j d9 d9' p(0) V{9 - 9') p{6') + O Q , (6.11) 

and hence the relationship (4.17) in the limit (3. 4). 7 

7 In deriving (4.17), we did not use the condition n — v/2 — > oo in (3.4), which amounts to requiring 
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Figure 1: The distribution p^\0) of the gauge field moduli in the PWMM around the background 
(3.3) is plotted for 7=1,2, 3, 4, 5, 16 with k — 16 and n = v = 31 at a temperature corresponding 
to x = 0.104 near the critical point x c ~ 0.072. (A similar behavior is obtained for decreasing 
I = 31, 30, • • • , 16.) The statistical errors are omitted since they are smaller than the symbol size. 
The solid line represents the result (2.10) for the J\f = 4 SYM on R x S 3 at the same temperature. 

7. Monte Carlo simulation of the effective theory 

In this section we confirm the large- N reduction by performing Monte Carlo simulation of 
the effective theory (4.3) for the gauge field moduli with the particular background (3.3). 
For the simulation details including the adopted algorithm, we refer the reader to ref. [48], 
where the effective theory (4.3) for simpler backgrounds was studied. In fig. 1 we plot the 
distribution p( T \0) of the gauge field moduli near the critical temperature for various /. 
As one goes towards the midpoint I = {y + l)/2, the distribution converges rapidly to the 
result (2.10) for the J\f = 4 SYM on R x S 3 at the same temperature. 

Thus we have seen that the distribution of the gauge field moduli has the property 
described below (6.9). From the argument in the previous section, this implies that the 
free energy of the PWMM agrees with that of SYM as in eq. (4.17). In order to confirm 
the agreement more explicitly, we calculate the free energy by using (4.16), where S e g is 
now replaced by the expectation value of the effective action (4.3) obtained by the Monte 
Carlo simulation 

We have performed simulations for k = 40, 100, 200 and n = v = 15, 23, 31 at x = 
0.071,0.074,0.077, which are near the critical point x c ~ 0.072 . We first extrapolate our 
results to k = oo. As is explained in ref. [48], the leading finite k effect for the normalized 
free energy is given by ^fp, which comes from the first term of (4.8). In fig. 2 (left) we 
plot our Monte Carlo results for the normalized free energy a 

gainst _|_ for n = v = 31. 

One can see that our data can be nicely fitted by a straight line. This allows us to make an 

the smallest fuzzy sphere to be regarded as a continuous sphere. It could therefore be that the equivalence 
holds irrespective of how one takes the large-n and large-i/ limits as far as n — > is satisfied. 
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extrapolation to k = oo. In fig. 2 (right) we plot the extrapolated values against 1/n, and 
find that they lie on a straight line. This allows us to extrapolate our results to n = oo. 
Thus we obtain the normalized free energy in the limit (3.4), which is plotted in fig. 3. Our 
results obtained from the PWMM agree nicely with the known result (2.12) for the Af = 4 
SYM on R x S 3 . 




Figure 2: (Left) The normalized free energy is plotted against -^f— for n = v = 31. The error bars 
represent the statistical error. The data can be nicely fitted to a straight line, which enables us to 
make a large-fc extrapolation. (Right) The results after large-fc extrapolation are plotted against 
1/n. The error bars represent the fitting error associated with the large-fc extrapolation. The data 
can be nicely fitted to a straight line, which enables us to make a large-n extrapolation. 




Figure 3: The normalized free energy of the PWMM around the background (3.3) is plotted 
against the dimensionless parameter x representing temperature near the critical point x c = 0.072. 
The data points are obtained by extrapolating results for k = 40, 100,200 and n = v = 15,23,31. 
The error bars represent the fitting error associated with the large-n extrapolation. The solid line 
represents the result (2.12) for the TV = 4 SYM on R x S 3 . 



8. Generalization to R x S 3 /Z q 



Our calculation can be generalized to M = 4 SYM on R x S 3 /Z q (q / 1) [51] and = 8 
SYM on R x S 2 [52], from which one obtains the PWMM by shrinking S 3 /Z q and S 2 
to a point, respectively. Unlike the M = 4 SYM on R x S" 3 , these SYM theories have 
only SU(2|4) supersymmetry. They have many classical vacua, all of which preserve this 
symmetry. The gravity dual corresponding to each vacuum is proposed in ref. [51]. The 
large-iV reduction is expected to work for each vacuum by choosing the corresponding 
background of the PWMM as one can see from refs. [32, 33, 34]. This regularization 
preserves the full SU(2|4) supersymmetry of the original theory for each vacuum. In what 
follows we present the calculation for M = 4 SYM on R x S 3 /Z q , which reduces to what 
we have seen above by setting q = 1. A simpler case of = 8 SYM on R x S 2 is discussed 
in appendix B. In both cases we restrict ourselves to the trivial vacuum for simplicity. 
Let us consider the theory (3.1) around the background (3.2) in the case 

kj = k, ni = n + q[l— UJ [ ) for I = 1, ■ ■ ■ , v , (8.1) 



and take the large- A limit in such a way that 8 

qv 

n ^ oo, v — > oo, fc ^ oo, n > oo 

2 

a 2 k 

with A = - — fixed . (8.2) 
n 

The resulting theory is expected to be equivalent to the planar limit of N = 4 SYM on 
R x S 3 /Z q [32, 33, 34] with the radius of S 3 and the 't Hooft coupling constant given, 
respectively, by 

2 

R s 3 = - > ^RxS A /Z q = ^ V S3/Z q , ( 8 -3) 

where V s z/z q = 2tt 2 (R S 3) 3 jq is the volume of S 3 /Z q . In what follows we confirm this 
statement in the weak coupling limit at finite temperature. The effective action is given 
by the same form as (4.10), where nj in is given by (8.1). 

Let us repeat the calculation in section 5 for general q. We first take the n — > oo limit 
in (8.2). Then, (4.13), (4.14) and (4.15) are reduced to 

zV> J \x) =x^- | ~- — ) , (8.4) 

4 J ' J) (z)=^-(- t— .:;--) - TV.: ! : t ■:>--!■ ^ 



(8.6) 




8 Strictly speaking, one has to tune A as a function of the UV cutoff parameters n and v according to 
the coupling constant renormalization in the non-conformal case q 7^ 1. This issue is irrelevant, though, in 
the weak coupling limit we are discussing. 
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Since these are Toeplitz matrices, we can set z^' J \x) = zf J \x) (i = s, v, f). The Fourier 
transforms (5.6) of these functions are evaluated as 

x(l + x 2 )(l - x 2q ) 2qx c ' +1 (u-2x cl + ux 2 i) 

Zsi * X ' ' ~ (l-x 2 ) 2 (l + x 2 <? -2xiu) + (1 - x 2 )(l + x 2 i - 2x q u) 2 ' *• ' 

x 2 (l+x 2 ) 2x q (u-x q )(q-l-{q-3)x 2 )+x 2 (3-x 2 )(l-x 2q ) 
z ^ x ' A ^ - (1-^2)2 + (l-x 2 ) 2 (l+x 2 <?-2x<?u) 

2qx 2q (2u 2 - 1 - 2ux q + x 2q ) + 2qx q+2 (u - 2x q + ux 2q ) 
+ (1 - x 2 )(l + x 2q - 2x q u) 2 ' [ ' 

2s § 2x q+1 2(u - x q )(2x 2 + g(l - x 2 )) + 2x1(1 - x 2q ) 
Z/(X ' } ~ (1 - x 2 ) 2 + (1 - x 2 ) 2 (l + x 2 <? - 2x q u) 

2qx 2q+ ^(2u 2 - 1 - 2ux q + x 2q ) + 2qx q+ l {u - 2x q + ux 2q ) 
+ (l-x 2 )(l+x 2 9-2x<?u) 2 ' ( ' 

where u = cos A. One can easily see that 



x(l + x 2 )(l - x 2q ) + 2qx q+1 (l - x 2 ) 



maxi s (x,A) = z a (x,0) - (1 - x 2 ) 2 (l - x?) 2 

A/ x , A , 4x 2 (l-x 2g )-2x«(l-x 2 ) 2 (l-x 9 ) + 2^x' ? (l-x 4 ) 
max^(x, A) = z„(x,0) = (1 _ x 2 )2(1 _ x9) 2 ' 

A . x . „ , . 2x5(1 + x)(x(l-X 2 «) + ^(1 -X 2 )) 

max%(x,A) = %(x,0) = ^ (1 _ x 2 )2(1 _ J )2 ' (8 ' 10) 

where the right-hand sides coincide with the single-particle partition functions in SYM on 
R x S 3 /Z q obtained at one loop in ref. [53]. By applying the same argument as in appendix 
A, we find 9 

nn\pV p (X)\ = Vi(0) . (8.11) 



mm 

A 



Therefore the critical temperature is determined by V±(0) = 0, which coincides with the 
condition for M = 4 SYM on R x S 3 /Z g . Thus we find that the critical temperature in 
N = 4 SYM onfix S" 3 /Zg is reproduced from PWMM. 

We can evaluate the free energy with an assumption made in section 6 and see that 
it agrees with the continuum theory. We have also performed Monte Carlo simulations of 
the effective theory for the gauge field moduli as we did for the q = 1 case and confirmed 
the assumption as well as the agreement of the free energy explicitly. 

Let us see what happens if we omit fermions and consider a bosonic theory. If we 
introduce r adjoint scalars instead of 6, the effective action takes the same form as (4.10) 
except that (4.12) has to be replaced by 

V^'J) = l{6u - rzi I ' J \x p ) - 4 7,J V)} • (8-12) 



9 For general q, we had to assume that the critical point lies in the regime i C 1 to derive (8.11) 
analytically. 
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Correspondingly, (5.7) has to be replaced by 

oo 

V P W= E V p (K) e tKX = -{l-rz s (xP,X)-z v (xP,X)} , (8.13) 

K=-oo P 

where z s (x p ,X) and z v (x p ,X) are given by (8.7) and (8.8), respectively. From (8.10), we 
find that (8.11) holds also for (8.13). Therefore, we see that the critical temperature of the 
bosonic YM theory on R x S 3 /Z q is reproduced from the corresponding bosonic matrix 
model. We have also checked the agreement of the free energy. 

9. Summary and discussions 

In this paper we have provided a nontrivial test of the large- N reduction for N = 4 SYM 
on R x S 3 , which enables us to study the AdS/CFT correspondence from first principles. 
By expanding the PWMM around a background corresponding to R x S 3 , we reproduced 
the deconfinement phase transition in M = 4 SYM on R x S 3 in the planar limit at weak 
coupling. The planar M = 4 SYM on R x S 3 can thus be regularized by the PWMM 
around that background. This regularization preserves the SU(2|4) supersymmetry, which 
is expected to enhance to the full superconformal SU(2,2|4) symmetry in the large-TV 
limit. Considering that we are dealing with a finite temperature set-up, which breaks su- 
persymmetry, our results actually suggest that the SO(4) symmetry, 10 which is the bosonic 
subgroup of the superconformal symmetry, is restored in the large- N limit. 

By choosing a different classical vacuum of the same model, we also reproduced the 
deconfinement transition in SYM on other space-times such as RxS 3 /Z q and RxS 2 . These 
theories only have SU(2|4) supersymmetry, which is fully preserved by our regularization. 
The gravity duals are known [51], and if one can calculate various quantities on the gravity 
side, one can study the gauge-gravity duality from first principles. It is interesting that the 
PWMM can be used to regularize all these theories in the planar limit by just changing 
the background configuration. 

We were also able to reproduce bosonic theories on curved spaces from the PWMM 
without fermions around the corresponding backgrounds. However, we consider that su- 
persymmetry is needed to protect the classical backgrounds against radiative corrections at 
strong coupling. Therefore, we consider that one has to choose sufficiently low temperature 
or to compactify the Euclidean time direction super symmetrically in order for the large- N 
reduction to work at strong coupling. 

Monte Carlo simulation of the PWMM can be done in exactly the same way as in the 
case of DO-brane system, which simply corresponds to the \i = case of the PWMM. In 
particular, in order to respect supersymmetry maximally, we consider it important not to 
discretize the Euclidean time direction, but to use finite numbers of Fourier modes after an 
appropriate gauge fixing [6, 8]. Indeed in the case of DO-brane system, the gauge-gravity 
duality has been confirmed with high precision including a' corrections [10, 11]. 

Restoration of the SO(4) symmetry was also checked by explicit one-loop calculations at zero temper- 
ature [32]. 
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The novel large- A r reduction discussed in this paper enables us to extend these numer- 
ical works to J\f = 4 SYM on R x S 3 in a straightforward manner. One only has to deform 
the one-dimensional gauge theory by the mass parameter fj,, and prepare an appropriate 
initial configuration to obtain R x S* 3 in the large- N limit. We consider it remarkable that 
there exists a seemingly feasible way to simulate N = 4 SYM on R x S 3 at strong cou- 
pling, and hence to investigate the AdS/CFT correspondence from first principles. Note, 
in particular, that the existing checks of the AdS/CFT correspondence are restricted to 
quantities protected by supersymmetry somehow. We hope that our formulation enables us 
to calculate quantities on the gauge theory side without such restrictions and to compare 
the results against predictions from the gravity side. 
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A. Derivation of eq. (5.11) 

In order to derive eq. (5.11) in section 5, let us note first that Zi(x, A) in (5.7) can be written 
explicitly as 

3 

„ , 2x 2 (l + 2u 2 - 4xu + x 2 ) , k , 

* (g ' A)= \l + x 2 -2xu) 2 (A - 2) 

where u = cos A. For instance, (A.l) can be derived as 

*.(*, A) =* ^ (* ( 1 + £ e dn*-HA)* + g e ^-* )K \ \ {AA) 



d ( x 



(A.5) 



dx \ 1 + x 2 — 2x cos A , 
Then one can easily find that 

max Zi(x, A) = zi(x, 0) = Zi(x) for i = s,v,f, (A. 6) 

A 

where Zi{x) is defined for the N = 4 SYM in eq. (2.7). For odd p, this implies that 



mm 

A 



m[ P V p (\)} = pV p (0) = pV p , (A.7) 
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where V p is given by eq. (2.5). Taking the derivatives of Zi(x), we find 



d 1 + 4x + x 2 

fa Zs{x) = (1-x) 4 



d . , 12x 

9 , , Q(x^ + xf ) /& \ 

= > • (A - 8) 

which imply pV^, > Vi for odd p / 1. Therefore, 

min {pF p (A)} = y 1 (0) = Vi . (A.9) 

A, p:odd I- ) 

The remaining task is to show that 

pVp(X) > Vi for even p . (A.10) 

Let us note first that 

d * -24x p {x p - x 3p ) + 8x§ p (l - x p )(l + Ax p + x 2p ) - 8x 2p {l - x p ){l - x^fu 
" (1 + x 2p - 2x p u) 3 ' 

(A.ll) 

from which we find that pV p (X) is minimized either at u = 1 or at u = — 1 for each p. Next, 
from (A. 8) we find that 

fi T P _i_ 1 2r 2p — 2r 3p 4- 1 firl p 
V 1 < 1 - 6z s (x p ) - z v (x p ) - 4z f (x p ) = 1 - + LZX {1 _ £ + = F p , (A.12) 

where we have defined a new function F p of x. Considering that 

pV p (0) -F p = 8z f (x p ) > , 

« , , „ 16x p (3x p + 2x 2p + 3x 3p + x^ p (l + x p ) 3 ) /K , 

- F P = ^ (1 _ x2p)3 ^ ^ > , (A.13) 

we obtain the inequality (A. 10). From (A.9) and (A. 10), we obtain (5.11). 
B. Large- iV reduction for Af = 8 SYM on R x S 2 

In this appendix we consider N = 8 SYM on R x S 2 . The "large- A reduction" in this case 
is nothing but the well-known construction of planar field theories using fuzzy spheres. We 
discuss it here nevertheless to see how our calculations reduce in this simpler case. 
Let us consider the theory (3.1) around the background (3.2) in the case 

v = 1, k\ = k, n\ = n , (B.l) 



- 18 - 



and take the large- N limit in such a way that (See footnote 8.) 



g 2 k 

n — > oo, fc — > oo with A = fixed . (B.2) 

n 

The resulting theory is equivalent to the planar limit of N = 8 SYM on R x S 12 , with the 
radius of S 2 and the 't Hooft coupling constant given, respectively, by 11 

R S 2 = - , A i?x5 2 = XV S 2 , (B.3) 
A* 

where V52 = 47r(i? >5 2) 2 is the volume of S" 2 . After taking the limit (B.2), we find that 
(4.13), (4.14) and (4.15) are reduced to 

4-r 2 

^ 1,1)(X) = (T^F ' (B ' 5) 

.d,i)^_ 2x 3 / 2 (l + x ) 
(1-x 2 ) 2 



(*) = 71 ; • (B-6) 



where the dimensionless parameter x is defined in eq. (2.6) with R S 3 = 2//x. Rewriting 
(B.4)~(B.6) in terms of 

x = exp(-^)=exp(-f)=* 2 , (B.7) 

they completely agree with the single-particle partition functions in J\f = 8 SYM on R x S 2 
obtained at one loop in ref. [50]. Therefore, the free energy agrees with that of M = 8 
SYM on R x S 2 . 

We can redo the calculation omitting fermions in the PWMM. The free energy of the 
resulting bosonic matrix model around the background (B.l) agrees in the limit (B.2) with 
the corresponding bosonic theory on R x S 2 . 
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